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ABSTRACT 

We study radiative gravitino decay within the framework of R- violating supersymmetry. 
For trilinear R-violating couplings that involve the third generation of fermions, or for light 
gravitinos, we find that the radiative loop-decay G ^ 'yu dominates over the tree-level 
ones for a wide set of parameters. We calculate the gravitino decay width and study its 
implications for cosmology and collider physics. Slow-decaying gravitinos are good dark 
matter candidates, for a range of parameters that would also predict observable R-violating 
signatures in colliders. In general the branching ratios are very dependent on the relative 
hierarchies of R-violating operators, and may provide relevant information on the flavour 
structure of the underlying fundamental theory. 
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1 Introduction 



The possibility of R- violating supersymmetry [1,2] has been extensively studied, as an 
alternative scenario to the Minimal Supersymmetric Standard Model (MSSM) [3]. Indeed, 
the symmetries of the Standard Model, in addition to the Yukawa couplings that generate 
fermion masses, allow for additional dimension-four couplings of the form 

XLiL^jEk + X'LiQ.Dk + X'ViDjDk (1.1) 

where the L{Q) are the left-handed lepton (quark) doublet superfields, and the E {D, U) 
are the corresponding left-handed singlet fields. Amongst these so-called R-violating cou- 
plings, the first two violate lepton number, while the third violates baryon number, and 
the symmetries of the theory imply that there are 45 operators in total — 9, 27, and 9 
respectively for the three terms, due to SU (2) and SU{'i) invariance. 

The stricter bounds on R-violating operators come from proton stability. However, it 
is by now well-understood that R-parity [4] is not the only symmetry that can guarantee 
proton stability; baryon or lepton parities [5, 6] can have the same effect. In fact, it is 
crucial to exclude the simultaneous presence of only certain products of LQD and U DD 
couplings [7]. In addition to proton stability, most operators are subject to experimental 
constraints from the non-observation of modifications to Standard Model processes, or of 
possible exotic processes [sjll. 

Going beyond proton stability and parameter constraints, one of the reasons that R- 
violating supersymmetry has not seemed very appealing, is the notion that in these models 
the lightest supersymmetric particle (LSP) is not stable, and therefore one has to search 
elsewhere for dark matter. 

However, how absolute is this statement? Is it possible that the lightest supersymmetric 
particle decays so slowly that cosmologically it is almost stable, while at the same time 
there is at least one coupling large enough to lead to observable collider signatures for 
R-violation? In this letter we shall take this to mean that at least one coupling is larger 
than approximately 10~^ for sparticles of 100 GeV mass. This could for instance happen, 
if the LSP is a light gravitino whose R-parity violating decays to conventional particles are 
for some reason very suppressed. Scenarios with light gravitinos are well-motivated — and 
one could even expect masses as low as 10~^ eV [10]. 

Tree-level R-violating gravitino decays have been studied in [11], where gravitinos were 
found to be able to decay before the present epoch for large values of the R-violating 
couplings and gravitino masses, but not earlier than the start of big-bang nucleosynthe- 
sis (BBN). It was concluded that the considered scenario would upset the standard BBN 
scenario, and did not seem to constitute a natural solution for the cosmological gravitino 
problem. Two-body gravitino decays to photon and neutrino, for bilinear R-parity vi- 
olation, were studied in [12] — with the decays arising from the mixing of neutralinos 
with neutrinos — and it was concluded that gravitinos of less than 1 GeV mass can be 

^Additional constraints can be obtained from a possible detection of photon or lepton plus missing 
energy signatures from NLSP decays to a gravitino LSP [9] . 
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the dark matter and at the same time consistent with neutrino mass generation from bi- 
hnear R- violation. Recently, in [13], these decays were revisited, and in models where 
R-parity breaking is tied to B-L breaking, very small R-violating couplings were predicted, 
in conjunction with a photon flux that can explain the apparent EGRET excess in the 
extragalactic diffuse gamma-ray flux [14] for ~ 10 GeV gravitinos. In addition, radiative 
neutralino decays in R-violating schemes have been studied in [2,15]. 
Here we will focus on radiative gravitino decays 



in schemes with explicitly broken R-violation from trilinear terms, and will compare them 
to the tree-level three-body decays that are expected from the same terms. Which processes 
will dominate, and what will be the cosmological and collider signatures depends on the 
flavour structure of the R-violating operators, which we will also discuss. 

In principle one may imagine that gravitinos are almost-stable and could "act" as dark 
matter in the following cases: 

• Suppose the only relevant operators are the \"U^DjDk. In the three-body gravitino 
decays discussed in [11] there would inevitably be a top-quark in the final state. If 
the gravitino is lighter than the top quark — even if heavier than all other fermions 
— then it is stable with respect to tree-body decays, up to mixing effect^. 

• This brings us to the next step: Suppose that all dominant operators involve third- 
generation fermions, and lead to bottom-quark or tau final states. For a gravitino 
below the ~ 1 GeV scale, tree-level decays will again be very suppressed. 

• Finally, for "super- light" gravitinos, as in [10], gravitinos are essentially stable under 
the 3-body decays! 

In either of the above cases — classified according to the possible range of gravitino 
masses for the suppression of three-body decays — gravitinos have a very large lifetime, 
that can exceed the age of the universe. In the case of a X'U^DjDk operator, tree-level 
decays can proceed only via quark mixing (which in some models could be very suppressed). 
For LLE or LQD operators, the radiative decay of the gravitino can still dominate over the 
tree-body modes. In particular for the cases where third-generation fermions are involved, 
where both the loop factor of the radiative mode, and the phase-space suppression of the 
tree- level diagrams are larger. 

In what follows, we will calculate the amplitude for radiative gravitino decays and dis- 
cuss the effects of the underlying flavour structure in Section [2|, leading to a comparison 
between radiative and tree-level gravitino decays in Section [31 We outline the related cos- 
mological implications for dark matter, photon spectra and big-bang nucleosynthesis, and 
discuss possible collider phenomenology from the decays of the next-to-lightest supersym- 
metric particle (NLSP) in Section H] that could help to further distinguish between the 
different scenarios of slowly decaying gravitinos, before we Conclude. 

^Note that for an operator of the form X'LiQ^Dk this argument does not hold, since, when we pass 
from superfields to component fields the LiQ^ part can become iit or Vib. 
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2 Radiative gravitino decays 



There are three main structures of diagrams that generate the decay f ll.2p . as shown in 
Fig. [H Note that the class of operators that can be involved are of the form LLjEj or 
LQjDj. In the latter case the neutrino is coupled to down-type quark to preserve SU(2) 
invariance. The common j index appears since the same fermion flavour should couple to 
both the gravitino and the photon. In Fig. [2] we show the corresponding three-body decay 
induced by the same operators. 



(a) 7 




k 



Figure 1: Basic set of Feynman diagrams for radiative gravitino decay, shown for (s)lepton 
loops. Arrows denote flow of lepton number for left-chiral fields. 




Figure 2: Three-body decay of gravitino via R-parity violating coupling. 
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The i?-parity violating coupling for incoming left-chiral fields can be parametrised a^ 



k'ui : -iXPL, (2.1) 

with A a dimensionless quantity, and Pl/r = |(1=F75)- The photonic coupling to sparticles 
is given by 

i+tA" : +ie{q+ - (2.2) 
with and the incoming slepton momenta. The gravitino couplings are (see [16,17]): 

CmA" : ^^Pfi7'^7^ (2.3) 

with M = (87rG'Ar)"^/2 = 2.4 x 10^^ GeV the reduced Planck mass. 

In order to have a non-zero contribution, there must be an even number of 7 matrices 
between the projection operators at the vertices where the gravitino and the neutrino 
couple. Thus, we immediately see that the amplitude is proportional to the lepton mass, 
rrii — or a down-type quark mass, for the LQD case. This amounts to a mass term 
insertion, as is required in order to get the correct helicity combinations at each vertex. 
In addition to the diagrams shown in Fig. [1] there is a corresponding set of diagrams with 
reversed fermion flow and right-handed (s)leptons. 

2.1 Feynman amplitudes and gauge invariance 

Extracting one lepton mass factor, the total amplitude for the Feynman diagrams shown 
in Fig. [1] can be written as 

Using the equations of motion of the gravitino [16]: 

'y^r{p) = 0, p,r{p) = 0, {^-m^)r{p) = 0, (2.5) 

and of the neutrino: ^u{p~k) = ^u{p — k) (neglecting the neutrino mass), one can simplify 
the matrix elements considerably. The dimensionless quantity J-' can thus be written as 

= u{p-k)[ak^^+b k^{e ■ p) + ce^] ^^{p). (2.6) 

The constraint of gauge invariance imposes a relation on a, h and c, and we are left with 
two gauge-independent structures. This is different from the case of xo 1^ [2,15] where 
there is only one amplitude. 



^For outgoing left-chiral fields, this becomes +iXPfi. 
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Using the kinematical relation, (p — k)"^ = 0, or 2p ■ k = m^, it is convenient to choose 
the two amphtudes as follows: 

J' = uip- k)Pn [ci^i) + c^A'j^^rip), (2.7) 

where Ci and C2 are dimensionless coefficients, the overall factor Pr arises from the coupling 
to the left-handed neutrino, and the structures A^f^^ and A^^^ are given by 

4^) = A;^[2(e-p)-/m^]/m|, 

A^ = {2k,{e-p)-e,ml]/m%. (2.8) 

The integrals from diagrams (a) and (b) will be logarithmically divergent and effectively 
regularised by the contributions of diagram (c). Thus, we see that can be expressed in 
terms of two three-point functions, and finite differences of two-point functions. These 
terms must separately be gauge independent: 

J^ = u{p- k) [XaCt'^ + X,Cj') + r(i)A5j') + F(2)A5(2)]^M(p)^ (2.9) 

where the coefficients X and Y must be of the forms given in (12. 8p . and the Cq are 
three-point functions corresponding to diagrams (a) and (b), whereas the A5o are finite 
differences of two-point functions. In the notation of LoopTools [18, 19], we have 

C'f^ = Co(m|, 0, 0, m|, m^, mj), 

Cq ^ = Coim^, 0, 0, m|, m'j, m|), 
Ai?^^^ = 25o(m|, m|, mj) - Bo{0, m|, rnj) - -Bo(0, mj, mj), 
A5f = 5o(m|, m|, mj) - 5o(0, m|, mj). (2.10) 

Adding the contributions of the three diagrams, we find that the coefficients of the 
amplitude are given by: 

Ci = 2[(m| - + m^)CQ"^ + 2i?o(m|, mj, mf) - Bo{0,mj,m'j) - _Bo(0, m^, m^)], 

C2 = 2[m2ci"^ + m| C^'^ + 5o(m|, m|, m^) - Bo{0, m|, mj)]. (2.11) 

These are both UV- finite, and have the structures given in eqs. (12. 9p and (l2.1Up . 
2.2 Flavour considerations 

Before comparing the radiative and the three-body decays, we would like to comment 
on the dependence of our results on flavour physics. The relative magnitude of radiative 
gravitino decays as compared to the tree-level ones depends on the flavour structure of the 
R-violating operators — as we shall see, for higher generations the radiative decay becomes 
larger and the same is true for the suppression of the tree-level diagrams. This picture is 
not unknown to us: the Yukawa couplings that generate fermion masses also have clear 
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hierarchies, and progressively lead to higher masses as we pass from the first to the third 
generation of quarks and leptons. 

Given that the same fermion fields that enter into R-violating operators, also enter 
in the Yukawa mass terms, one may try to directly link R-violating hierarchies to those 
of fermion masses [20,21]. This is done using models with family symmetries. Fermion 
generation charges are chosen in such a way that only the third generation mass terms 
have a zero charge (and thus are allowed when the symmetry is exact). The rest of the 
masses are generated at a higher order by the spontaneous breaking of this symmetry by 
the vacuum expectation values of singlet fields, and are suppressed by the heavy mass scales 
of the theory [22,23]. If R parity is violated in such models, couplings with different family 
charges will also appear with different powers of the family symmetry-breaking parameter, 
and thus with different magnitudes. 

In general, for models appearing in the literature, the relative flavour charges and thus 
the mass matrices, are determined by the GUT multiplet structure; particles in the same 
GUT multiplet have the same charge. Nevertheless, in all cases, the observed fermion hier- 
archies require smaller charges for the operators of the higher generations — typically zero 
for the top Yukawa mass terms, but also for the bottom and tau in supersymmetric models 
with large tan/5. This implies that it is reasonable to expect that R-violating operators 
that contain fields of the third generation may also be larger. Whether this happens will 
depend on whether fermion hierarchies can be directly linked to the R-violating ones, or 
whether extra singlet fields with non-zero flavour charge are involved, thus reversing the 
hierarchies with respect to those of the masses [21]. In this latter case, tree-level gravitino 
decays are more likely to dominate, unless gravitinos are very light. 

Moreover, one should appropriately take into account mixing effects. Indeed, even in 
the case of one dominant operator in the basis of current eigenstates, mixing effects will 
induce non-zero coefficients for related operators in the basis of mass eigenstates. These 
will be suppressed by the mixing parameters with respect to the dominant operator, but 
will be non-zero, and this may affect phenomenological and cosmological prediction^. We 
should also keep in mind that experiments only provide information on the Cabibbo- 
Kobayashi-Maskawa (CKM) quark mixing matrix V'"^^ = V^'^V^, and that one can 
construct theoretical models where the left quark mixing is either in the up sector, or in 
the down sector, or both. Similarly, the Maki-Nakagawa-Sakata lepton mixing comes from 
the product of those of charged leptons and neutrinos; with the additional complication 
that, for the latter, we have both Dirac and Majorana mass terms. 

Without entering into detailed model building, we generically observe the following: 

(i) The right-handed quark mixing (relevant for U and D) is essentially not constrained 
by the data. Therefore, for a model with left-right asymmetric mass matrices, one 
could also imagine a manifestation with minimal mixing in the right-handed sector, 

"'The fact that there are strict bounds on some operators, imphes that mixing effects may in given 
models generate additional bounds on couplings that at a first glance look less constrained. This has been 
analysed in detail in [21], where it was shown that in theories with strong correlations between operators 
(such as left-right symmetric models), the effects can be particularly significant. 
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in which case a U^DiDj operator would be the only relevant one, and the gravitino 
would be essentially stable! 

(ii) For the left quark mixing (relevant for Q), we know the values from Vckm, where 
for instance the 2-3 mixing is a factor of ~ 0.04. Thus, a coupling X'L^Q^D^ also 
implies the coupling 0.04A'L3Q2-D3. 

(iii) The left lepton mixing (relevant for L) is constrained by the lepton data, giving large 
1-2 and 2-3 mixing, and small 1-3. 

We see that there are in principle several flavour choices that can lead to significant 
effects from the decays under discussion, and one could in fact reverse the argument: the 
study of the relative magnitude of radiative versus tree-level violating decays, may yield 
relevant information for the flavour structure of the underlying theory. 



3 Radiative versus three-body decays 

The decay rate for the radiative decay G ^ 7Z/ is given by 

1 1 



IGvr niQ 



\M\\ (3.1) 



where \ denotes the absolute square of the Feynman amplitude, averaged over gravitino 
spin states, and summed over photon polarisations. In terms of the decomposition (12.41) . 
the decay rate can be given by 



G I -Tr|2 



20 48774 M2 

with 



^1, (3.2) 



4 

spin, pol. 

^|ci|' + ^|c2|'-^Re(ctc2), (3.3) 



where the averaging over gravitino spin states involves 
PAp) = ^'^^l{p)^u{p) 

spin 

-i*- { - ^) - i - ^) - ^) V} . (3,4) 

The diagrams with reversed fermion flow, and right-handed s(fermions) in the loop are 
found to have the same structure, differing only in the sfermion loop mass. The decay rate 
for G ^ 7z/ is identical, and will provide a factor of two in the total rate. 



8 



In Fig. [3l we compare the resulting lifetime of the gravitino for radiative versus tree-level 
decays resulting from trilinear R-violating terms. We fix = 200 GeV and A = 0.001. 
The interested reader should rescale the results to suit his or her preferred value of A, 
using r oc 1/A^. The black lines denote radiative decays for (s)tau and (s)muon loops, 
while the green lines represent the results for the corresponding three-body decays, with 
intermediate stau and smuon, taken from [11]. 




mg [GeV] 

Figure 3: Gravitino lifetime versus mass for two-body (black) and three-body decays 
(green). Also shown is the area excluded by a universe age of 13.7 Gyr (blue), and the 
kinematical threshold for the two-tau final state of the three-body decay (red line). 

We see that the radiative decays can easily dominate for gravitino masses below 30 GeV. 
While the three-body decay involving an intermediate stau hits the kinematical threshold 
at 2mr, the radiative dominance is still present for decays involving an intermediate smuon, 
where there is no threshold at the mass scales shown. For a selectron intermediary (not 
shown) the conclusions are the same, with radiative dominance below a gravitino mass of 
around 2 GeV. This behaviour is controlled by the mass dependence of the decay width: 
for the three-body decay F^. oc m^, while for the radiative decay Tq oc rriQ at low gravitino 
masses. This can be understood from Eq. (12.111) . where we see that the leading term in ci 
and C2, for the limit of slepton masses much larger than the gravitino and lepton masses, 
2mjCo, is essentially independent of mass. The physical interpretation of this is that the 
gravitational coupling compensates for the high loop-mass by its increasing strength for 
higher masses. Because of the helicity structure of the couplings, the two-body decay 
width is also heavily dependent on the fermion mass, oc mf at low gravitino masses. Thus 
dominant third generation couplings, which imply a tau or bottom quark, give significantly 
shorter lifetimes for the radiative decay. 
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To constitute a realistic dark matter candidate the gravitino lifetime should exceed the 
age of the universe. This lower boundary is shown as a blue area in Fig. [31 For the R- 
parity violating couplings considered here this only allows us to exclude gravitino masses 
above ~ 150 GeV. However, we shall see in the next Section that the photon flux from the 
gravitino decays can provide us with a much stronger bound. 

In Fig. m we show the dependence of the gravitino lifetime on the stau mass for the 
sum of radiative and three-body decays, and again for a fixed A233 = 0.001. As expected 
the gravitino mass where the three-body decay takes over depends on the stau mass. The 
insensitivity of the lifetime to the slepton mass for light gravitinos can again be understood 
from the dominant terms in Ci and C2. 




Figure 4: Gravitino lifetime versus mass for different slepton masses. 



4 Cosmological implications and collider signatures 

From the above we see that the mass and lifetime of the gravitino in principle allow it to be 
dark matter. Whether this could be the case for a specific model depends on the predicted 
gravitino density. If the universe has gone through a period of inflation, the primordial 
gravitino abundance is erased, and gravitinos are created either thermally or from the 
decays of the NLSF^. This implies that the abundance of gravitinos is very sensitive to the 
reheating temperature of the universe — which is an open question — and can well be in 
the correct range for gravitinos to be dark matter, or a component of it. 

''With non-zero R-parity violation the NLSP production mechanism no longer applies. 
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In MSSM dark matter considerations, the allowed parameter space is severely con- 
strained by possible effects of late NLSP decays on the light element abundances, as pre- 
dicted by Big Bang Nucleosynthesis [24]. However, as has been remarked by several au- 
thors [12,13,25], in models where the NLSP decays via R- violating couplings, the BBN 
bounds are easily satisfied. 

With regards to baryogenesis and leptogenesis it has been observed that severe con- 
straints on i?-parity violating interactions can be derived from cosmological arguments [26] : 
a pre-existing baryon asymmetry would be erased if the interactions are strong enough to 
come to equilibrium at the time of the electroweak phase transition. However, it has 
been argued that this need not be the case in various schemes where flavour effects in 
leptogenesis/baryogenesis are considered [27]. In our analysis therefore we consider larger 
R- violating couplings than those used in e.g. [13]. 

Our final cosmological consideration is whether the presence of photons in the final 
state is compatible with measurements on the extragalactic diffuse photon background, 
and may even be reconciled with an apparent excess in the EGRET data, as proposed 
in [13]. This will in fact give the strictest bounds on R- violating couplings, for the case of 
radiative decays. 

The photon flux from gravitino decays is described by a red-shifted monochromatic 
line, which is given by [13]: 

FiE) = = BR(G ^ 7z/) x CJl + Kx^)-^/^x'>/^e(l - x), (4.1) 

dE 

where 

X = , = -T- and k = ——. (4.2) 

With current values for the cosmological parameters [28] we arrive ai|^ 

/ in^^s \ 

= 1.06 • 10"^ GeV cm^2 sr^^ s"^ and /t ~ 3. (4.3) 

For comparison, the original EGRET analysis [29] gave a power law description of the 
extragalactic flux as 



.dJ /I GeV 



0.1 



E''^ = 1.37 ■ 10~^ ( ) GeV cm'^ sr^^ s-\ (4.4) 

dE \ E J 

in the energy range 30 MeV to 100 GeV. The non-observation of an excess with monochro- 
matic origin can be used to restrict the gravitino mass and lifetime. 

In Fig. [5l we show the maximum allowed value of the R- violating coupling A for a range 
of gravitino masse^. We require that the gravitinos can be dark matter, with a lifetime of at 



^Note that there is an order of magnitude difference between our value, and that of [13]. We have learnt 
from the authors that this is due to a misprint in their paper. 
^Similar considerations are directly applicable to A'. 
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least 10 times the current age of the universe, and that their radiative decays are consistent 
with the photon spectrum measured by EGRET, extrapolating the power law behaviour 
up to 400 GeV. The result is a simple linear dependence on the log of the gravitino mass 
for the range where the radiative decays dominate. For higher gravitino masses, where the 
three-body decays become important, the bound relaxes as the branching ratio to photons 
decreases dramatically. Eventually the bound on the lifetime, for gravitinos to be dark 
matter, becomes dominant and the maximum allowed value decreases rapidly towards the 
slepton mass threshold. The extra features seen at the very highest masses are due to 
threshold effects for the three-body decays. 




Figure 5: Maximum value of R-violating coupling A versus gravitino mass. 

We find that the allowed couplings are small for radiative decays in an intermediate 
mass range, but can be significantly larger when tree-level decays dominate or for very light 
gravitino masses. In all cases the couplings can be sufficiently large to lead to interesting 
expectations for collider phenomenology. 

It is interesting to note that the strict upper bounds on the couplings derived from 
the photon spectrum, naturally bring us to the range that is relevant for the generation of 
neutrino masses [2, 8, 30] compatible with neutrino data. Moreover, the power law excess 
claimed for the 2-10 GeV photon energy range in the more recent analysis of the EGRET 
data [14] could easily be explained by decaying gravitino dark matter with a mass of around 
10 GeV and R-violating couplings of C(10~'^). 

Moving on to the collider signatures, these will be the standard signatures of R-violating 
super symmetry, with multi-lepton or multi-jet events and the possibility of explicit lepton 
number violation at the final state. Given the magnitude of R-violating couplings required 
from the photon spectrum and for gravitino dark matter, MSSM production of sparticle 
pairs followed by R-violating decays would be expected. However, for the parameter space 
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and flavour structure where tree-body decays dominate, single superparticle productions 
may still occur along the lines previously proposed in the literature [1,2,31]. 

The specific decay modes depend on the NLSP and the branching ratios of R- violating 
versus R-conserving modes. If the NLSP is the neutralino it will decay to three fermions, 
with diagrams similar to those of the tree-level decay of the gravitino, see Fig. O If instead 
the NLSP is a slepton — e.g. z/ or even — its decay will dominantly proceed through 
a direct decay to two fermions, for the cases with large R-violating couplings involving 
its fiavour, or through an intermediary neutralino to four fermions (for small R-violating 
couplings). The specific branching ratios depend on the supersymmetric parameter space 
and the magnitude of the R-violating coupling, and have been extensively studied [1,2,31, 
32]. 



NLSP 


LLE 


LQD 
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X? 
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Table 1: NSLP R-violating decays. 



The choice of NLSP in the above table has been made according to the principle that, for 
universal sparticle masses at a high scale, the running of couplings to low energies implies 
an increase in masses due to gauge interactions and a decrease due to Yukawa couplings. 
Thus, the neutralino and the sneutrino, having no electric charge, and the right-handed 
stau, due to the lack of weak interactions and a large Yukawa coupling in supersymmetric 
models with large tan (3, are the most obvious candidates. 

5 Conclusions 

We have studied radiative gravitino decays to a photon and a neutrino, generated from 
trilinear R-violating operators. We calculated the decay rate, and compared it to tree-level 
decays of the gravitino to three fermions, which arise from the same operators. There 
is a wide region of parameters for which the loop suppression of the radiative mode is 
less effective than the phase-space suppression of the tree-level one, particularly for light 
gravitino masses and R-violating operators that involve the third generation of fermions. 

Whether enhanced gravitino decays can be expected, depends on the flavour structure of 
the R-violating operators. Without entering into details of model building, it is nevertheless 
possible to make simple analogies, and link R-violating couplings to those that generate 
Yukawa mass terms for the fermions. Mixing effects can also be important, since tree-level 
decays that are disallowed by phase-space considerations for certain operator flavours could 
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nevertheless proceed through quark mixing effects, suppressed by the mixing factors. For 
a dominant operator U^DiDj, the absence of the radiative mode and the top quark in the 
vertex, imply that the gravitino can be essentially stable in models with small right-handed 
quark mixing. 

We find that the gravitino lifetime is typically longer than the age of the universe due to 
suppression by Planck scale and loop efi^ects, meaning that gravitinos can constitute dark 
matter. The R- violating couplings and sparticle masses in such scenarios have values that 
can give rise to distinct signals in colliders, and the R-violating couplings are in general 
higher than those expected in bilinear R-violating schemes [13]. 

However, for the parameter space where radiative decays dominate, photon spectra 
further constrain the magnitude of R-violating couplings, particularly for and 
Li 2fiQ?,D?, operators. In fact the presence of photons in the final state of R-violating 
gravitino decays is compatible with an excess seen in the EGRET data on cxtragalactic 
diffuse gamma-rays, as proposed in [13]. Moreover, the range of couplings allowed by these 
considerations is compatible with the one required in order to generate acceptable neutrino 
mass patterns. 

Finally, the decays of the next-to-lightcst super symmetric particle can still proceed via 
the R-violating vertexes at significant rates, and are thus easily accommodated within the 
framework of Big Bang Nucleosynthesis. 
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